Pumping through a quantum dot in the proximity of a superconductor 
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We study adiabatic pumping through a quantum dot tunnel-coupled to one normal and one 
superconducting lead. We generalize a formula which relates the pumped charge through a quantum 
dot with Coulomb interaction to the instantaneous local Green's function of the dot, to systems 
containing a superconducting lead. First, we apply this formula to the case of a non-interacting, 
single-level quantum dot in different temperature regimes and for different parameter choices, and 
we compare the results with the case of a system comprising only normal leads. Then we study 
the infinite-U Anderson model with a superconducting lead at zero temperature, and we discuss the 
effect of the proximity of the superconductor on the pumped charge. 

PACS numbers: 73.23.-b, 74.45. +c 



I. INTRODUCTION 

A finite charge can be pumped through a meso- 
scopic system in the absence of an applied bias volt- 
age by changing periodically in time some parame- 
ters of the system. If the parameters are changed in 
time slowly with respect to the lifetime of the elec- 
trons in the system, pumping is adiabatic. The idea 
of electron pumping through a mesoscopic system is 
based on a work of Thoulessi Since then a large num- 
ber of thfioretjcal 2 i 3 i 4 i B i a i 7 i 8 i Q i 10 i 11 i ia i 13 t u i 16 i ia i ir i 18 i 19 and 
experimenta l 20 ' 21 ' 22 ' 23 works has been dedicated to this 
field. Adiabatic pumping has a geometric nature: The 
pumped charge depends only on the area spanned by 
the pumping cycle in parameter space and not on its 
detailed timing. It can therefore be related to geo- 
metric phases.— iSbLIS In normal conductors transport is 
due to the transfer of individual quasiparticles. In su- 
perconducting devices, quasiparticle transport is hin- 
dered by the gap, and charge pumping is due to Cooper 
pairs i 10 ' 11 ' 12 ! 13 In hybrid normal-superconducting sys- 
tems, Andreev reflection at the interface between the 
normal and the superconducting region allows for sub- 
gap charge transport in the system. In the absence of 
Coulomb interaction, adiabatic pumping in such hybrid 
systems can be studied by generalizing in Nambu space 
the scattering approach to adiabatic pumping i 14 i 15 ' 16 i 17 

In the present paper we study a system consisting of 
a quantum-dot attached to a normal conducting lead 
and a superconducting lead (N-dot-S system) as shown 
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FIG. I: Sketch of the dot attached to two leads. We are 
interested in the N-dot-S system, where the left lead is nor- 
mal conducting and the right lead is superconducting and will 
compare to the N-dot-N system, where both leads are normal 
conducting. Left and right lead are attached to the dot via 
tunnel barriers. The tunneling strength is described by the 
respective rate Fn and Fs. Possible time-dependent parame- 
ters are: the strength of the tunnel couplings Tn and Ts, and 
the dot-level position e. 

in Fig. [T] In the absence of Coulomb repulsion in the 
dot, adiabatic pum ping t hrough such a system has been 
discussed in Refs. HUT! On the other hand, when 
Coulomb repulsion is present, only DC transport through 
the N-dot-S system has been investigated! 24 ' 25 ' 26 ' 27 ' 28 ' 29 

We generalize a recently developed formalism to cal- 
culate the pumped charge through an interacting quan- 
tum do t 18 ' 19 to the case of the N-dot-S system. First, 
we apply this formalism to a non-interacting quantum 
dot, studying systematically the ratio between the charge 
pumped through the N-dot-S system and the charge 
pumped through a system consisting of a dot attached to 
two normal conducting leads (N-dot-N system). This ra- 
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tio was calculated in Refs. HUEi for a non-interacting 
resonant system at zero temperature. We consider finite 
temperature and extend the analysis to the off-resonant 
situation. Finally, we apply the formalism to the infinitc- 
U Anderson model in the Kondo regime. We discuss the 
pumped charge as a function of the level position and 
the superconducting gap and compare to the interacting 
N-dot-N system. 



II. MODEL AND FORMALISM 

In this paper we consider a quantum dot tunnel cou- 
pled to one normal and one superconducting lead, as 
shown in Fig. [T] The system is described by the Hamil- 
tonian: H = ^ a=NS (ff + H tumha ) + H dot . The dot 

Hamiltonian reads -ffdot = So- e (t)d\,d a + Un^n^, where 
a is a spin label, n a = d\d a and we allow the dot level 
e to depend on time t. The Hamiltonian for the normal 
lead reads = J2k a 6kC N ka c N ka' anc ^ the one f° r ^ ne 
superconducting lead reads Hg = ^2 k a CkC s ka c s ka + 



perconducting order parameter. We take A to be k- 
independent, which is appropriate for conventional BCS 
superconductors such as aluminum and niobium. The 
operator d\, (d a ) creates (annihilates) a dot electron 
with spin a, the operator ka (c a ^kcr) creates (anni- 
hilates) an electron in lead a with momentum k and 
spin a. The tunnel coupling between the dot and the 
leads is taken into account by the tunneling Hamiltonians 

H tU nn, a = Y,k,a ( V <* (*) c a,fco-^ + #- c -)> where a = N,S, 
and the tunneling amplitudes V a can be time-dependent. 
We assume that the phase of V a does not depend on 
time, since a time-dependent phase would describe an 
applied bias voltage. From now on we work in Nambu 
space, defining the fields: * Q , fe = (c a M , _ ki ), and 

$=(d v d\). 



E fe (Ac| 



S,feT c S,- 



fcj 



H.c), with A £ R being the su- 



Generalizing the approach of Refs. 18.19, the pumped 
current can be related to the local Green's function of 
the dot. In particular, the first adiabatic correction to 
the charge current from the dot to the left (normal) con- 
ductor through the barrier is given by 



'w~sM-U 



Tr 



N (*)G r ( W ,t)] (Gl(u,t)) X G8(a;,t) 



(1) 
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where f(u>) is the Fermi function, the caret indicates a 
matrix in Nambu space and T3 is the Pauli matrix given 

by t 3 = f J ^1 ) r ^^ 1C ^ rS ^ ^ Crm ™ ® corre- 
sponds to the so called average-time approximation, 18 
while J corr (discussed below) contains contributions due 
to vertex corrections^ and it is zero for all cases studied 
in this paper. The instantaneous retarded Green's func- 
tion is defined as G^(r,0,i) = -ie(r)({$+(r), $(0)}>, 
and its Fourier transform is Gq(uj, t) — J dr e LU1T G} ) {T, t). 
The last argument (t) of the Green's function is the 
time with respect to which the adiabatic expansion has 
been performed.— In particular, the instantaneous re- 
tarded Green's function, Go(w,t) is computed with all 
time-dependent parameters frozen at time t. The ad- 
vanced Green's function is related to the retarded one 
by Gg(w,i) = (Gg(w,i)) t . The intrinsic line width Y a is 
defined as T a (t, t') = 2np a V a {t)V*{t') with a = N, S and 
r Q (i) = T a (t,t). The normal lead is supposed to have 
flat bands with a constant density of states /On- The den- 
sity of states ps is the one of the superconducting lead in 
its normal state and it is constant, too. 

We now discuss the average-time approximation and 
give an expression for the correction term J c °" . The 
average time approximation neglects terms due to ver- 
tex corrections 19 in evaluating the adiabatic expan- 



sion of the self energy. The self energy is defined 
by the Dyson equation G(t',t",t) = g(t',t",t) + 
fdtifdt 2 G(t',t 1 ,t)Y,(t 1 ,t2,t)g(t 2 ,t",t) and it is a 
functional of the time-dependent Hamiltonian H(t), 
where r £ [t±, t%] . The adiabatic expansion is performed 
by linearizing the time-dependence of the Hamiltonian 
around a fixed time t (denoted as a third argument of 
the self energy). In order to solve the Dyson equation 
for the adiabatic expansion of the Green's function, the 
average-time approximation is performed on top of the 
adiabatic expansion. Such approximation consists in re- 
placing the linear time dependence on r by the average 
time (ti+t2)/2 of the interval [£1,^2] and thus it neglects 
the term 



E corr (t 1 ,t 2 ,t) = 

Tk f dr — 

Jk SH(t) 



H(r)=H{t) 



where the time integral is calculated along the Keldysh 
contour, Tk is the time-ordering operator along the 
Keldysh contour, and 5T,/5H (r) a functional derivative. 
Finally, the correction to the pumped current can be 
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written as 

jcovr 



(2) 



idbJ 



doj' 



ReTr 



cj' — cj — «0+ 

It is important to point out that the average-time 
approximation is exact whenever: i) the dot is non- 
interacting, ii) temperature is zero, the interaction is ar- 
bitrary, and the system can be effectively mapped to a 
Fermi liquid, and iii) when interaction is infinite and the 
self-energy £ is calculated up to linear order in T. There- 
fore the correction term equals zero in all cases discussed 
in this paper. 

In the following we are interested in calculating the 
charge pumped through the system in the weak pumping 
regime (bilinear response in the pumping fields 30 ). The 
pumped charge Q is related to the pumped current by 
Q = J Q J (t), where T is the period of the cycle. We can 
choose as pumping parameters, X(t) and Y(t), any two of 
the three quantities e (t) = e + Se (t), T s (t) = f s +ST s _(t) 
and Fn (i) = Tn + ST-^ (t) . The time-averaged part X is 
denoted by a bar and the time-dependent part by 8X (t) . 
For the pumped charge, with pumping parameters indi- 
cated in the subscript, we find 



Q 



X,Y 



ev(X,Y) 



4tt 



iLj: I 

OOJ 
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^((fN-f s )^)^((^)- 1 Gg)}.(3) 



The prefactor 77 (X, Y) accounts for the amplitude and 
the relative phase of the pumping parameters and 
is equal to the surface in parameter space enclosed 
in one pumping cycle. It is defined as r\ (X, Y) — 

J dr [■§pSX(r)\ 5Y(t). The order of the parameters in 
the argument of r\ is important and has to be respected 
in the formulas above (changing their order results in 
an additional minus sign.) Gq and Gg are the retarded 
and advanced dot Green's functions (both matrices in 
Nambu space although the caret has been omitted to 
simplify the notation) with all parameters taken at their 
time-averaged value. 



III. NONINTERACTION G QUANTUM DOT 

We start by considering the simple case of a non- 
interacting quantum dot. The results which we obtain 
for the non- interacting dot, can equivalcntly be calcu- 
lated by a generalization of Brouwer's formula to the N- 
dot-S systemi 14 i 15 i 16 The instantaneous Green's function 



of the dot is 



(4) 



where the self energy due to the proximity of the super- 
conducting lead reads 



Es = 



VA 2 -u 
A 



r s 



2 VA^Z 



■ (5) 

\ 2 VA 2 -^ 2 VA 2 -^ / 

We now focus on the limit of large superconducting gap 
A -> 00, when (S s )i,2 = (Ss)a,i = r s /2 and (S s )i,i = 

(S S )2,2 = 0. 

Let us consider first the particle-particle spectral func- 
tion, which is defined as Ai.i(uj) — — (l/ir)Im{G\ A. 
The spectral function conveys important information on 
the proximity effect in the quantum dot. Furthermore, 
the DC transport properties are influenced by the energy 
distribution of spectral weight, for example the Andreev- 
reflection probability R^(oj) has a similar structure as 
Ai,i(u>). Due to the close relation of the spectral func- 
tion with the retarded Green's function appearing in the 
formula for the pumped charge, the spectral function will 
also prove to be a useful quantity to analyze pumping in 
the N-dot-S system. It turns out that Ai t i is given by a 
linear combination of two Lorentzians and reads 



Ai,i(w) 



+ 




,(6) 



with the Lorentzian bell defined by Lr(uJ — E) = 

(57)/ (i w ~ ^") 2 + (^) 2 ) • P resence 01 two peaks in 
the particle-particle spectral function [Eq. ((6])] is a signa- 
ture of the proximity effect induced by the superconduc- 
tor. The width of the Lorentzians appearing in Eq. ([6]) is 
equal to the coupling to the normal lead, while their po- 
sitions are shifted, with respect to the dot-level position, 
by the coupling to the superconducting lead Ts- This 
fact can be referred to as a proximization of the dot-level 
position, which might be contrasted to the level renor- 
malization due to Coulomb interaction. The repulsion of 
the two Lorentzian peaks with increasing Tg reflects the 
increasing coupling between electron and hole excitations 
in the dot [represented by the off-diagonal terms in Eq. 
©]■ 

In the non-interacting case Eq. ([3]) can be written in 
a compact form if the two pumping parameters, denoted 
by X (t) and Y (t), are chosen among the three quantities 
{If - (t) , (t) , e (t)}M. We denote by Z the parameter 
which is kept constant in time. Eq. ([3]) simplifies to 



4e 

}x,y = —r](X,Y)e x ,Y,z 

7T 



(7) 
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FIG. 2: Ratio between the pumped charge in the N-dot-S and 
in the N-dot-N systems, Qx,y IQ^x^ > as a function of f s/Fn- 
The curve is the same for any couple of pumping parameters 
X, Y . Temperature is equal to zero and the bare dot level is 
at resonance. 



FIG. 3: Ratio between the pumped charge in N-dot-S and 
in the N-dot-N system, Qx y/Q^x Yi as a function of Ts/Fn- 
Temperature is equal to ksT = 1.5Fn. The pumping param- 
eters are Tn and Ts (solid line), Tn and e (dashed line) and 
Ts and e (dotted and dashed line). The bare dot level is at 
resonance. 



where N is the denominator of the particle-particle spec- 
tral function with time-averaged parameters 



N = 



2 



(8) 
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The quantity r\ (X, Y) is the surface in parameter space 
enclosed in one pumping cycle. The antisymmetric ten- 
sor Ex,y,z has to be evaluated respecting the order of 
the parameters (t) , ^# (t) , e (i)}. In the case of 
weak-pumping with the two tunnel-barrier strengths, i.e 
X(t), Y{t) = f s,n + £Ts,n W, an d Z — e, we can rewrite 
Eq. as 



^ N ,r s = ^(r N ,r s ) i i-|^(-g 



df \ ORa 
de • 



(9) 



which relates the pumped charge to the probability of An- 
dreev reflection Ra, computed with the time-dependent 
parameters frozen at their average value. The probability 
of Andreev reflection is given by^ 



Ra 



F 2 F 2 
1 N 1 S 




+ rg,(r§ 



4e 2 ) 



r- 
1 s 



(10) 



We are now interested in comparing our results to 
the ones obtained for an N-dot-N system (see Fig. [1] 
with normal-conducting right lead) . The current pumped 



through a non-interacting dot attached to two nor- 
mal conducting leads can be obtained using Brouwer's 
formula^ For the case where the tunneling rates Fn(£) 
and Fs(t) are the pumping parameters, it has been found 
that the current in the presence of one superconducting 
lead is about a factor 4 bigger than in the case of two nor- 
mal leads, if the dot level is at resonance (e = 0) and tem- 



perature is zero 
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We find exactly the same behavior 



for all three choices of pumping parameters, at zero tem- 
perature. It turns out that the ratio QxyIQxy ( ec l ua l 
for any pairs of pumping parameters) has a maximum, as 
a function of the level position, at resonance. 33 We plot 
in Fig. rj] the ratio Qxy/Qxyi f° r e = 0, as a function 
of Ts/T]$: a maximum is present when Tq/Tjsi — 0.9. 
In the limit of small and large Ts/Tn the charge ratio 
goes to zero. This happens because, in the presence of a 
superconducting lead, pumping is mediated by the An- 
dreev reflection probability, Eq. (fTTJ|) . On one hand, for 
Ts/Tn *C 1, the ratio is suppressed since Ra implies 
higher-order tunneling processes to the superconducting 
lead with respect to tunneling to the normal lead. On the 
other hand, for Ts/Fn 3> 1, the Lorentzian contributions 
to Ra get narrower with decreasing Tn and more distant 
with increasing Ts, thus suppressing the Andreev reflec- 
tion probability at lu — (relevant for zero temperature). 
In the case of finite temperature the charge ratio depends 
on the parameter choice and can take values much bigger 
than 4. Fig. [3] shows the ratio Qx y/Qx y f° r the three 
different parameter choices at temperature k^T = 1.5Fn, 
with the dot level being at resonance. 

In the following we discuss some qualitative features, 
in the high-temperature limit (k^T 3> Tn), of the ratio 
Qx y/Qx y- Expressions for the pumped charge are re- 
ported in Appendix [Al While for the N-dot-S system the 
leading order of the pumped charge is the first deriva- 
tive of the Fermi function, in the N-dot-N system this is 
the case only when the level position is one of the pump- 
ing parameters. The first contribution to Q 



N 

r N ,r s 



starts 
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with the second derivative of the Fermi function. On 
the other hand, the parameter Tg enters in the spectral 
density of the N-dot-S system via the proximized peak 
positions. The high-temperature limit of the charge ra- 
tio as a function of temperature is therefore a constant if 
{r Q , e} are the pumping parameters and it is quadratic in 
temperature in the case that {Tn^s} are the pumping 
parameters: 



Ts 



Q N 



(k B TY 



TnI^s 



rN,rs 



f 2 r 2 \ • 

4 4 J 



i f N f g (f N + f s ) ; 

r 2 r 2 

^4^4 



(11a) 



(lib) 



rs, 



1 f N f s (f 



Q N 




For all three cases described by Eqs. (fTT|) . the ratio of 
the pumped charge, as a function of the average dot level 
position, has a maximum when the average dot level is 
at resonance. At e = 0, with Fg = Fn, we find that 
QkJQr N ,e - 2 and Qr^/QfL - 4. Furthermore 
we note that due to the second term in Eq. (jllcp . the 
ratio Qp s e /Qr s c is increasing in a monotonic way for 



increasing r s /r N , while <2r N ,r s /Qr N ,r s and Qr> 



have a maximum as a function of Ts/Tn- This behavior 
is present in Fig. [3J where we also observe a large value 
of the maximum of Qp N rg /Qr N r s > which is due to the 
quadratic dependence on k^T . 

The plots in Fig. [4]show the pumped charge, as a func- 
tion of the average dot level position, through the N-dot- 
S system [panel (a)] and through the N-dot-N system 
[panel (b)] for different choices of pumping-parameter 
pairs at k^T = Fn and for Fn = Fg. We point out sev- 
eral differences: In the N-dot-N case the pumped charge 
is the same for pumping with {rN,e} and with {Fs,e}. 
In the presence of a superconducting lead the symmetry 
between the two leads is absent as soon as temperature 
differs from zero. Furthermore we report that the width 
of the pumped charge is strongly temperature dependent 
for the N-dot-N system, while for the N-dot-S system the 
width is saturating for high temperatures to a value that 
depends solely on the coupling to the leads, Tn and Tg. 
This is due to the fact that the gap of the superconductor 
is always bigger than temperature and no quasiparticle 
transport, which is sensitive to temperature, takes place 
between the superconductor and the dot. The symmetry 
of the pumped charge as a function of the dot-level posi- 
tion is present for both systems. For the N-dot-S system, 
when pumping with Tn and Tg , this happens despite the 
fact that, at least at high temperatures, pumping is dom- 
inated by the proximized level position. 
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FIG. 4: Pumped charge in the N-dot-S [panel (a)] and in the 
N-dot-N [panel (b)] systems as a function of the average level 
position e. Temperature is set to k^T = Fn and Fn = Fs. 
The pumping parameters are {Fn, Ts} (solid line), {Fn, e} 
(dashed line) and {Fs, e} (dashed-dotted line). The units on 
the x-axes are the ones determining the width of the respective 
function. 



It might be interesting to contrast the pumped charge 
in the N-dot-S system with the DC linear conduc- 
tance of the same system (here all parameters are time- 
independent): 



4e 2 

G = ir 



(12) 



In the high temperature limit, i.e. k^T 3> Tn, the con- 
ductance is given by 



G 



4e 2 7T r N r 



h 4£| 



Nl S 




(13) 



IV. STRONG INTERACTION LIMIT 

In the regime of strong Coulomb repulsion double oc- 
cupation is forbidden and coherent Andreev scattering 
can occur only if A is finite (electrons can enter the su- 
perconductor within a time window of order h/ A). In the 
Kondo limit, when temperature is zero, we can evaluate 
the retarded equilibrium dot Green's function making use 
of the mean-field slave-boson technique ! 26 ' 34 As the dot 
can only be empty or singly occupied we can describe 
the singly-occupied dot state by the fermion operator f a 
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FIG. 5: Pumped charge Qr N ,r s [panel (a)] and Q%v s [panel 
(b)] as a function of the average level position for differ- 
ent values of A , solid lines, and Qr N ,r s and Qr N , E i dashed 
lines. Temperature is zero and we consider symmetric barriers 
(f N = f s ). The cut-off energy is W = 20f N . 



(known as pscudo fcrmion) and the empty dot by the bo- 
son operator b. The real dot operator d a is simply given 
by d a = w f a . The operators b and f a need to fulfill the 
constraint b^b + J^a fa fa = 1- Such a constraint is intro- 
duced in the Hamiltonian using the Lagrange multiplier 
A which plays the role of a chemical potential. In the 
mean- field limit, the boson operator b is replaced by a 
real number bo- The quantities bo and A are determined 
by minimizing the free energy. Note that we are allowed 
to use an equilibrium formalism, since after performing 
the adiabatic approximation we are left with a current 
formula, Eq. fl}, that contains only equilibrium quanti- 
ties. 

Details of the derivation of the retarded dot Green's 
function by means of the slave-boson technique, in the 
limit A > Tg, are given in Appendix [Bl It is important 
to notice that the pumped charge in the deep Kondo 
regime is zero. Therefore, in the following, we will be in- 
terested in the mixed-valence regime only. The pumped 
charge in the weak-pumping limit for the aforementioned 
three pairs of pumping parameters can be calculated us- 
ing Eq. The pumped charge has a similar structure 
in all three cases. We find, in particular, the general 
relation 



1 



1 
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(14) 




which is also valid for the noninteracting N-dot-S system 
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FIG. 6: The peak position e pca k is shown as a function of 
the superconducting gap for different pairs of pumping pa- 
rameters on a logarithmic scale. The pumping parameters 
are {r^Ts} (solid line) and {FN,e} (dashed line). The tem- 
perature is zero and Fn = Ts = 0.05W. 



at zero temperature. 

In analyzing the properties of the pumped charge as a 
function of the different system parameters we can profit 
from the fact that, in the Kondo regime, the system be- 
haves like a non-interacting system with renormalized pa- 
rameters, denoted as Fn, Ts and e. We should however 
keep in mind that the renormalized parameters depend 
in a complex way on all the bare parameters. Panel (a) 
of Fig. shows the pumped charge Qr N r s ( s °hd lines) 
as a function of the level position, for increasing values of 
the superconducting gap A. The charge Qp N Fs in the N- 
dot-N is also shown for reference (dashed line). In panel 
(b), the pumped charge is shown for Tn and e chosen as 
pumping parameters. We note that for both the N-dot-N 
system and the N-dot-S system (for any value of the su- 
perconducting gap A) the charge has a bell-like structure. 
The width of the bell-like curves reflects the width of the 
respective particle-particle spectral functions, which for 
the N-dot-S system is given by Tn while for the N-dot-N 
system it is given by Tn + Ts . For the values of parame- 
ters considered in Fig. [5] (f n = f § = 0.05VF), the curve 
relative to A = 2Tn is still very far away from the curve 
relative to the N-dot-N system. 

The striking similarity of the pumped charge for the 
two different pairs of pumping parameters [compare Figs. 
[5] (a) and (b) and Eq. (|T3])] can be explained by the 
fact that pumping takes place via renormalized parame- 
ters. Indeed, in the zero-temperature regime with strong 
Coulomb interaction, the Green's function of the system 
depends on the bare parameters solely via renormalized 
parameters r^Ts and e. Therefore, the detailed behav- 
ior of the pumped charge is determined by the derivatives 
of the renormalized parameters with respect to the bare 
ones. In particular when the bare dot level is deep below 
the Fermi energy, i.e. in the Kondo regime, the renormal- 
ized parameters are independent of the bare parameters 
and therefore the pumped charge is zero. 

Now, we focus our attention on the dependence of the 
pumped charge on the superconducting gap. In a non- 
interacting dot, the pumped charge is always independent 



7 




0.1 0.2 0.3 0.4 0.5 



AAV 

FIG. 7: Peak width Ae pca k for the pumped charge plot- 
ted as function of the superconducting gap A. The pump- 
ing parameters are {Fr^Ts}. The temperature is zero and 
f N = f s = 0.05W. 

of the superconducting gap A at T — 0, since the Green's 
function of Eq. ([JJ, for uj = 0, is independent of A. The 
situation is different in the presence of Coulomb interac- 
tion, where the pumped charge depends on the supercon- 
ducting gap via the renormalized parameters. Now, we 
investigate how the features of the bell curves, describing 
the pumped charge as a function of the average level e, 
depend on the gap A. 

Peak position. Fig. [5] shows that the position of the 
peak of the pumped charge Q s depends on In A. The 
peak position is shifted towards zero (from below) for in- 
creasing gap. This is consistent with the fact that the 
renormalized level acquires a logarithmic correction in A 
proportional to 

Peak width. Fig. [7] shows the peak width Ae pca k (half- 
height-width) for the pumped charge Qp N Fs as a func- 
tion of the superconducting gap A. As stated before, the 
peak width has an analogous behavior to the width of the 
spectral function. The coupling to the superconducting 
lead influences the width of the spectral function via the 
renormalized coupling Tn- The decrease of Tn with in- 
creasing superconducting gap induces the A-dependence 
shown in Fig. [71 

V. CONCLUSIONS 

We have extended a Green's function formalism to 
compute the pumped current through a quantum dot 



with Coulomb interaction, to the case when the dot is 
tunnel-coupled to one normal and one superconducting 
lead. First, we have applied our formalism to a non- 
interacting dot in an N-dot-S configuration, where we de- 
rived a compact formula for the pumped charge, and we 
calculated the ratio between the charge pumped through 
the N-dot-S system and through an N-dot-N system. At 
finite temperature, we found different characteristic be- 
haviors of this ratio, depending on the choice of pump- 
ing parameters. Then, we have studied the pumped 
charge through a quantum dot with infinitely strong 
Coulomb interaction at zero temperature using the mean- 
field slave-boson approach. In particular, we have fo- 
cused our attention on the influence of the superconduct- 
ing gap A on the pumped charge. We have explained 
features of the pumped charge exploiting the mapping 
to a non-interacting system with renormalized parame- 
ters. In the end we wish to comment on the experimental 
realization of the pumping mechanism discussed theoret- 
ically in this article. The system depicted in Fig. [JJ can 
be realized in analogy to an experiment reported recently 
by van Dam et al£^. Here a semiconductor quantum dot 
is contacted with superconductors, using indium arsenide 
nanowires combined with local gate voltages to create a 
quantum dot with tunable coupling to the leads. Su- 
perconductivity is induced in the wires by means of the 
proximity effect using aluminum contacts. In order to 
realize pumping, time-dependent control of the system 
parameters is required. 



APPENDIX A: PUMPED CHARGE IN THE 
NON-INTERACTING DOT: k B T > r N 



In this appendix we report the expressions for the 
pumped charge through the non-interacting dot in the 
weak-pumping regime in the limit k^T 3> Tn- 
N-dot-S system From Eq. ([JJ, for the three possible 
choices of pumping fields, we obtain: 



?r N ,r s - er? (F N , r s ) p 



4 F 2 I r N I r s duj 

4 4 
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N-dot-N system We compare the results for the N-dot- 
S system to the charge pumped through the N-dot-N 
system by varying the respective pumping parameters. 
We find: 



Qr N , e 



Q N 



e . d 2 f 

2^(rN,r s )^(e) 



-2e?7(r N ,e) 



3/ 



(f N + f S ) 

r N 9/ 



r N + r & 



2 da; 



9a; 



(A4) 
(e) (A5) 

(A6) 



where the trace is performed in Nambu space and ui n = 
7r(2n + l)//3 is a fcrmionic Matsubara frequency. The 
pseudo-fermion Green's function Q has the same func- 
tional form of the noninteracting Matsubara Green's 
function Q, but contains the renormalized parameters 
Tn = ^o^n, Tg = b^Ts and e = e + A. Minimizing the 
free energy with respect to the variables A and bo we find 
the equations: 



oo 

= b% + k B T Tr{g(iov)a z } (B2a) 



APPENDIX B: MEAN-FIELD SLAVE-BOSON 
TECHNIQUE FOR THE N-DOT-S SYSTEM 

In this Appendix we report details of the slave-boson 
method for the N-dot-S system^ We calculate instan- 
taneous equilibrium Green's functions where all param- 
eters are always taken at their time-averaged value. We 
therefore omit the bar and, e.g., write Tn instead of Fn- 
Furthermore all Green's functions considered in this ap- 
pendix are 2x2 matrices in Nambu space and we omit 
the caret. 

The free energy of the system can be written as a func- 
tion of the Matsubara Green's function Q [related to the 
retarded Green's function by the analytic continuation 



oo 

F = -k B T J2 Trjlng- 1 ^)} 



Xbi 



= X + k B T J2 Tr{0(io; n )r(ia; n )} . (B2b) 

n— — oo 

The matrix in Nambu space T (iu> n ) is given by 



r (mj„) 



-i^sign (uj n ) 









-i^sign(o;„) 



2 J^l+^' 2 
r s a 



2 Ju>l+A 2 

a r s 



2 y^f+A 3 



(B3) 



Equations (|B2a ) and (b) are used to determine A and 
&o. At zero temperature we can replace the sum over the 
Matsubara frequencies by an integral. Furthermore, we 
consider the limit A ^> Ts and we perform the following 



e, (Bl) approximation 



2(1 



n— — oo 



r (iw„) 



-i^sign (u n ) 



2 



Es 

2 



-i-^sign (w p 



) sign (w n ) 





* + if) si s n K) 



for a; < A 



for uj > A . 



(B4) 
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Under such an approximation, Eqs. (|B2b ) and (b) read 



^5 



A = 




(B5a) 



(B5b) 



-- N-dot-S 




FIG. 8: Spectral density of the dot in the N-dot-N (solid 
line) and in the N-dot-S system (dashed line) in the Kondo 
regime (e = -0.25W, T N = T s = OMW). The gap of the 
superconductor is half as big as the bandwidth cut-off (A = 
0.5W). The Kondo-temperature of the N-dot-N and the N- 
dot-S system differ by a huge factor: T k N = 7.8 x 10" 4 W and 
T k s = f .3 x f(T 6 W. 



where we introduced a cut-off energy W representing the 
bandwidth of the leads. 

Finally, the retarded dot Green's function G rot is ob- 
tained from the retarded pseudo-fermion Green's func- 
tion G ret by means of the relation 

G rct = bl& e \ (B6) 

From the retarded Green's function of the N-dot-S 
system, we can calculate the spectral density Af 1 — 
— (l/ir)ImG r i\ and compare to the one for the N-dot- 
N system, in the same regime. The result is shown in 
Fig. [5] We find that the Kondo temperature in the 
N-dot-S system is much smaller than in the N-dot-N sys- 
tem (note that the width of the two curves is rescaled 
by the respective Kondo temperature). Nevertheless, in 
the deep Kondo regime, the spectral weight at the Fermi 
energy is the same for both systems. The strength of the 
Kondo-effect is therefore not affected by the presence of 
the superconductor in the limit of a large but finite gap 
A. 
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